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ABSTRACT 

Lyman-a spheres, i.e. regions around the first stars which are illuminated by Lyman a (hereafter, 
Ly a ) photons and show 21cm absorption feature against the CMB, are smoking guns at the dawn 
of the reionization epoch. Though overwhelming radio foreground makes their detections extremely 
difficult, we pointed out that, strong gravitational lensing can significantly improve their observational 
feasibility. Since Ly a spheres have ~ 10 sizes, comparable to the caustic size of galaxy clusters, 
individual images of each strongly lensed Ly a sphere often merge together and form single structures 
in the 21cm sky with irregular shapes. Using high-resolution ./V-body ACDM simulations, we found 
that the lensing probability to have magnification bigger than 10 is ~ 10 -5 . This results in > 10 6 
strongly lensed Ly a spheres across the sky, which should be the primary targets for first detections 
of Ly a spheres. Although the required total radio array collecting area for their detection is large 
(~ 100 km 2 ), the design of long fixed cylindrical reflectors can significantly reduce the total cost of 
such array to the level of the square kilometer array (SKA) and makes the detection of these very 
first objects feasible. 

Subject headings: cosmology: galaxy clusters - gravitational lensing 



1. INTRODUCTION 

In the concordance cosmology, first stars and galax- 
ies formed at redshift z ~ 30 in dark matter halos of 
mass > 10 6 Mq. Inefficient gas cooling, due to the ab- 
sence of heavy elements, possibly results in first stars of 
mass around several hundred solar mass. These first stars 
emit numerous near Ly a photons directly. Furthermore, 
secondary Ly a photons can be generated through the 
impact excitations by electrons produced by the X-rays 
from the high temperature stellar photosphere. Through 
frequent Ly a scattering in surrounding neutral hydrogen 
atoms, the neutral hydrogen spin temperature is coupled 
to its kinetic temperature. Since the kinetic temperature 
is lower than the CMB temperature at z ~ 30, a region of 
21cm absorption ag ainst the CMB develops around each 
first stars /galaxies (|Cenl l2006t IChen fc Miralda-Escude! 
2006). Such regions, named "Ly a spheres" or 21cm 
absorption "halos", are about of 10 -100 , much larger 
than the first stars/galaxies themselves and are possi- 
bly the first distinctive structures in the universe ac- 
cessible to observations. As unique features of first 
stars/galaxies, they will allow direct probe of the very 
first stage of the reionization process. This approach 
is quite complementary to the statistical approach of 
the 21cm observaton, such as measuring 21cm intensity 
power spectrum and bispectrum. It will shed light di- 
rectly on the reionization mechanism and could break 
degeneracies arising in the modeling of statistical quanti- 
ties, e.g., degeneracy in various reionization mechanisms. 

However, detections of these first structures are ex- 
tremely difficult, due to overwhelming radio foregrounds 
and stringent resolution requirement. For 21cm absorp- 
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tion "halos" a round stars in the 10 8 M Q dark matter ha- 
los (|CeiJl2006t l. a total collecting area of A to t ~ 100 km 2 
is required to resolve t hem and a At ot of ~ 10 km 2 is re- 
quired to detect them (|Zhang et al.l l2006). Ly a spheres 
around the first stars w hich formed in Kf'M^ dark mat - 
ter halos, as discussed in lChen fc Miralda-Escu de (2006), 

are even smaller (~ 10 ) with weaker signals. Their de- 
tection requires A to t <i 10 3 km 2 . Such mission would be 
extremely difficult. 

Fortunately, as shown in this paper, strong gravita- 
tional lensing of these Ly a spheres makes the first de- 
tection significantly easier. Galaxy clusters, the so called 
gravitational telescopes, have the power to magnify these 
objects by a factor of 10 or more and thus can sig- 
nificantly improve the observational feasibility. As es- 
timated in this paper, a survey capable of detecting 
these strongly lensed Ly a spheres can be built with 
cost comparable to the square kilometer array (SKA 4 ). 
Such survey will map the 21cm sky to very high accu- 
racy over a large fraction of the sky and a wide red- 
shift range. It will have significantly impact on both 
the fields of cosmology and reionization. For the cosmol- 
ogy part, it will allow precisio n lensing reconstruction 
through the 21cm backgr ound (|Cooravil200l iPenl 12004 
IZahn fc Zaldarriaga 20Qa) an d precision measurement o f 
the matter power spectrum (jLoeb fc Zaldarriagall2004j ). 
Furthermore, it will also resolve 21cm absorption 'ha- 
los', which again have powerfu l cosmological applications 
(jCenll2006l : IZhang et al.ll200l . 

Using strong l ensing to detec t high z objects is not 
a new idea (see iFort fc Mellierl (|1994 ) for an early re- 
view). Using this method, several hi gh redshift galax- 
ies at z ~ 6-10 ha ve been detected (jPello et "all 12004 
iRichard et alj I2006D . Ly a spheres lie at higher red- 
shifts, so the probability to be highly magnified is slightly 
higher. More importantly, strong lensing of Ly a spheres 
carry a feature distinctive to that of ordinary galaxies, 

4 http://www.skatelescope.org/ 
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explained as follows. 

Ordinary galaxies at cosmological distance typically 
have ~ 1 size, this size is usually much smaller than 
the caustic size of lens. So images of lensed galaxies 
are often separated. However, Ly a spheres have size 
~ 10 , such a size is comparable to the caustic size of 
galaxies clusters. As a consequence, individual images of 
each strongly lensed Ly a sphere often merge together 
and form a single structure in the 21cm sky with irregu- 
lar shapes. Although individual images of larger source 
usually has lower magnification, the merged image can 
still have a larger magnification. This effect reduces the 
requirement on A iot . Since there are hundreds of bil- 
lions of Ly a spheres, even a tiny probability to have a 
large magnification would result in a significant number 
of strongly lensed merged images. Interferometer array 
with Atot <C 1000 km 2 would be sufficient to detect these 
lensed Ly a spheres. 

The outline of this paper is as follows. In §2, we discuss 
the high resolution N-body numerical simulations we use 
and the method we employ to calculate the lensing prob- 
ability. In §3, we describe the strong lensing features of 
Ly a spheres and present the result of the lensing proba- 
bility. In §4 we estimate the observational requirements 
and feasibility and we summarize our results in §5. 

2. NUMERICAL SIMULATION AND LENSING METHOD 

The cosmological model considered here is the concor- 
dance ACDM model with the dimensionless matter den- 
sity fi m .o = 0.3, the cosmological constant J7a.o = 0.7, 
the shape parameter T = O m ,o/i = 0.21, er 8 = 0.9, where 
h is the Hubble constant in units of 100 km s" 1 Mpc -1 
and we take h = 0.7. A cosmological N-body simula- 
tion with a box size L = 300/i _1 Mpc, whic h was gener- 
ated with our vect orized-parallel P 3 M code (|Jing fc Sutol 
20021 Uingl [20021. is used in this paper. The simu- 
lation uses 512 3 particles, so the particle mass m p is 
1.67 x lO^fe^MfT y The gravitational force is softened 
with the S2 form (|Hocknev fe Eastwoodlll98ll ) with the 
softening parameter r\ taken to be iQhr 1 kpc. Since 
strongly magnified images are produced mostly at radii 
> 100/i _1 kpc in the lens plane, the resolutions arc suffi- 
cient (see also lDalal eh alJl2004 1. 

Dark matter halos are identified with the friends- 
of-friends method using a linking length equal to 
0.2 times the mean particle separation. The halo 
mass M is defined as the virial mass enclosed within 
the vir ial radius according to the spherical collapse 
model (iKitayama fe Sutol Il996t iBrvan fe Norman! 119981: 
Uing fc Sutol 12002ft . 

For a given cluster, we calculate th e smoothed sur - 
face density maps using the method of iLi et all (|2006h . 
Specifically, for any line of sight, we obtain the surface 
density on a 1024 x 1024 grid covering a square of (co- 
moving) side length of 6/i _1 Mpc centered on each clus- 
ter. The projection depth is chosen to be 6/i _1 Mpc. 
Notice that the size of the region is chosen such that 
particles within a few virial radii are included. Parti- 
cles outside this cube and large-scale structures do not 
contribute significantly to the lens ing cross-section (e.g. 
ILi et al.ll2005t iHennawi et al.ll2007h . Our projection and 
smoothing method uses a smoothed particle hydrody- 
namics (SPH) kernel to distribute the particle mass on a 



three-dimension grid and then integrat e along the line of 
sight to obtain the surface density (see lLi et al.ll2006l for 
detail). In this work, the number of neighbors used in 
the SPH smoothing kernel is fixed to be 32. We calculate 
the surface density along three perpendicular directions 
for each cluster. Once a surface density map is obtained, 
we compute the cross-section of highly magni fied images 
following the method given in ILi et all (f2005h . 

We put the sources at z s = 30, a typical red- 
shift of Ly a spheres. The dependence of strong 
lensing probability on z s at z s > 10 is negligi- 
ble, due to the weak dependence of source distance 
on z s and the fact that no clusters at redshift be- 
yond several. We simplify the background sources as 
spheres with diameter d s = 5", 7.5", 10", 15", 20". This 
choice of source size roughly c overed the source size in 
IChen fc Miralda-Escudi (|2006[) . Throughout this paper, 
we have neglected the difference in the 21cm brightness 
temperature across Ly a spheres. Since gravitational 
lensing preserves the 21cm brightness temperature, from 
the mapping between source position and lensing po- 
sition, the temperature distribution can be recovered 
straightforwardly. 

We generate a large number of background sources 
within a rectangle box. This rectangle is chosen to 
enclose all the high- magnification regions (/i > 1.7 and 
/i < 0) that can potentially form highly magnified im- 
ages. Such regions are much larger than the size of 
caustics and provide a reliable lensing cross section of 
image with /i > 2. The sources are located on a reg- 
ular grid that cover this rectangle with a resolution of 
0.5". For each source, ray-tracing is used to find the 
resulting image(s). The magnification of each image is 
defined as the ratio of the image and source's areas. The 
resolution in image plane is kept as 2" for sources with 
d s = 5", 7.5", 10" and 4" for d s = 15", 20", so the to- 
tal number of pixel is large than 40 for the image with 
magnification larger than 2. We calculate the total cross 
sections of the top 1000~2000 most massive clusters in 
each simulation output. Finally, we obtain the average 
cross section per unit comoving volume by: 



where (Ji(zi) is the average cross-section of the three pro- 
jections of the z-th cluster at redshift z\, and V is the co- 
moving volume of the simulation box. The optical depth 
can then be calculated as: 

where D s is the angular diameter distance to the source, 
z s is the source redshift, and dV p (z) is the proper volume 
of a spherical shell with redshift from z to z+dz. We used 
28 simulation outputs with redshift from 0.2 to 3. The 
integration step size is the same as the redshift interval 
of simulation output (dz « 0.1). 

3. STRONG LENSING PROBABILITY 

To illustrate the signatures of lensed Ly a spheres, we 
choose a cluster at z = 0.36 as the lens. By varying the 
size of the Ly a sphere, we find an interesting lensing fea- 
ture. Usually, the magnification decreases quickly with 
expanding source size, because the magnification which 
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Fig. 1. — The strong lensing signatures of first star Ly a spheres. 
Sources at z = 30 with different diameter d s = 2.5", 5", 10", 20" 
are shown in purple, blue, red and black respectively in the bot- 
tom panel. The corresponding images are shown in the top panel 
with same colors. The lens is at redshift z = 0.36 with mass, 
Mvir = 1-5 x 10 15 /i -1 Mq. The blue curves are critical curves 
and the green curves are caustics. When increasing the size of the 
source to cross the caustics, some separate lensed images will merge 
together. For a sufficiently large source, all lensed images will 
merge together and form a single structure with irregular shape. 
This feature is distinctive to ordinary strong lensing of galaxies. 
The resulting single lensed images will most likely be the first de- 
tected Ly a spheres. 



close to the caustic follows ~ \i oc y~ l , where y is the 
distance to caustic from inner. On the other hand, there 
are also two effects to increase the high magnification re- 
gion. When a source is close to but outside of the caustic, 
the magnification is modest for a point source, but for an 
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Fig. 2. — The optical depth as a function of magnification. Differ- 
ent color of line corresponding different size of source respectively. 
The optical depth for fi > 2 virtually does not depend on the source 
size. This is because the region with fj, > 2 in source plane is much 
larger than the size of caustic or the size of source. Hence, the 
cross-section with fj, > 2 is dominated by the low magnification re- 
gion and the optical depths are almost same for the point sources 
and the finite sources. 



extended source which is large enough to cover the caus- 
tic, a highly magnified image results. So the finite source 
can extend the high magnification region to the outside 
of the caustics. When a source is inside but faraway from 
the caustics, it will create multiple images with also mod- 
est magnifications, but if the source is large enough to 
match the caustic(s), some or all of these images will 
merge together, then the final magnification is still very 
high. Fig. [1] shows such merger effect of images. The lens 
is at redshift 0.36 with mass, M vir = 1.5 x 10 15 h~ 1 M e 
and the sources are fixed at redshift 30. Notice, there are 
no overlaps for a merged image. The incomplete images 
connect together along the critical curve(s) and build up 
such a larger image. 

Since the typical size of Ly a sphere is ~ 10 and 
is indeed comparable to the size of the caustics, it is 
frequent that lensed Ly a sphere will show up as a single 
merged image. The area of this image is a factor of fi of 
the source, where /z is the total magnification. 

Fig. [2] shows the optical depth as function of magni- 
fication. For typical size of Ly a spheres ~ 10 , there 
is 10~ 5 probability to produce /i > 10 merged images. 
Since there are hundreds of billions Ly a spheres, we ex- 
pect to find hundreds of thousands lensed Ly a spheres 
with /j, > 10. 

3.1. Model uncertainties 

We further investigate the completeness of our simu- 
lated cluster sample in mass and redshift range and thus 
the robustness of our predicted -P(^t). Lower panels of 
fig. [3] show the differential probability distribution for 
the optical depth as a function of lens redshift, for various 
source sizes and magnification thresholds. The optimal 
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lens redshift for larger source size and larger magnifica- 
tion threshold is lower. For the source size 5 ,10 , f 5 
and \x > 5, 10 investigated, the optimal lens redshift 
varies between 0.4 and 1. Thus cluster at redshift around 
0.4 to 1 will be most likely responsible for strong lensing 
of Ly a spheres. 

Top panels of fig. [3] show differential cross section 
as a function of logarithm lens mass, d<7/dlog(M)/<7 to tai) 
where M is the virial mass of individual lens cluster and 
"total is the summation of the mean cross-sections of all 
clusters. Since the optimal redshift is between 0.4 and 
1.0, we choose a median redshift z — 0.75. For this red- 
shift, we select the top 2000 most massive clusters in the 
simulation box. The low mass cut is about 3 x 10 13 and 
the corresponding number of particle within in virial ra- 
dius is ~ 2000. The small cluster has too small caustic(s) 
to highly magnify a big source. But when the source size 
is smaller, the source can have non-negligible probability 
to be strongly lensed by a less massive cluster, which may 
not be well resolved in our simulation. Meanwhile, when 
estimating the lens optical depth of small fj,, contribution 
from less massive clusters may be also non-negligible for 
their huge abundance. Thus, for small \i and small source 
size, our cluster sample can be incomplete. Top left panel 
of Fig. [3] shows that this is the indeed case for source size 
5 and [i > 5. However, for the more interesting case, 
10 and/or /i > 10, contribution from clusters less mas- 
sive than 10 13 ' 5 ft. _1 M Q is negligible. Thus, our cluster 
sample is complete at low-mass end for 10 sources. The 
last figure tells us that the optical depths are complete at 
low-mass end for /i > 10 for any size of source. Particu- 
larly, the lensing probability is 10~ 5 at this magnification 
threshold. 

On the other hand, Fig. [3] implies that clusters more 
massive than 10 15 ft~ 1 -M Q still have non-negligible con- 
tribution. Due to the finite box size, massive clusters 
are not fairly represented in our simulation. The lensing 
probability of /i > 10 may be under-estimated by a fac- 
tor no larger than 2 for the lack of very massive clusters 
in our simulation. 

As seen from Fig. [3J halos less massive than 
10 13 ' 5 /i _1 Mq have negligible contribution to ^ > 10 
lensing. Thus, we do not expect galaxies, which are not 
included in our simulation, can affect the above result 
significantly. This can be further understood as follows. 
For an isothermal distribution of matter in galaxies, the 
Einstein radius 6e = (cr„/186km s -1 ) 2 £>ls / D g arcsec. 
Since it is much smaller than the size of Ly a spheres, 
it is not surprising that they do not contribute much to 
strong lensing of Ly a spheres. 

We have to say that the lensing probability we give here 
is just a low limit in the concordance ACDM model. The 
baryon effect is not taken into account. In the real uni- 
verse, baryons account for roughly 20% of the total mass, 
which can cool (form stars) and sink toward the centers 
of clusters. The radiative cooling likely has two effects: 
it will increase the concentration of baryons at the center 
of clusters, and at the same time, make the clusters mor e 
spherical (e.g. iDubinskil 11994 [Kazantzidis et al.ll2004f ). 
The former increases while the latter decreases the lens- 
ing cross-sections, and so the overall influence dep ends 
on which effect d omina tes. iMeneghetti et al.l (|2003l ) and 
iPuchwein et alj (|2005t ) found that the increase in the 
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Fig. 3. — The two top panels show the differential cross section 
as a function of logarithm lens mass. The magnification thresh- 
old is labeled in each panel and the representing of color for each 
histogram is as same as in top and bottom panels. The total area 
under each curve is normalized to unity. The two bottom panels 
show the differential probability distribution for the optical depth 
as a function of lens redshift. The optimal lens redshift changes 
from 0.5 to 1 following the changing of source size and magnifi- 
cation threshold. Larger source size and magnification threshold 
causes lower optimal lensing redshift. 



lensing cross-section of giant arc due to baryon is quit 
modest, by a factor of < 2 for massive cluster and the 
the increase becomes larger for low-mass cluster. So, we 
can estimate the probability as ~ 2 x 10~ 5 if we take into 
account the baryon. 

From the two top panels in fig. [3j we can see our cluster 
sample is usually incomplete at low-mass end (for small 
source and low magnification threshold) or at high-mass 
end (for large source and high magnification threshold). 
This leads to an under-estimation of the lensing proba- 
bility by a factor of < 50%. In this paper, we use thin 
lens approximation but the large-scale structures along 
the line of sight also should have additional contribution 
to th e lensing probability although this is very modest 
(e-g 



h e lensing probability although this 
, iLi et al.|[2005t iHennawi et at| [2007) 



On the other hand, recent results from the three- 
year data of the Wilkinson Microwave Anisotropy 
Probe (WMAP) prefers lower cr 8 = 0.74 and Q mfi = 
0.238. than that in the concordance ACDM model 
(|Spergel et all l2006h . Since the abundance of massive 
clusters decreases dramatically in this low a% universe, 
the number of giant arcs is a factor of ~ 6 l ower than that 
in the concordance model (|Li et al. 2006b). Roughly, we 
expect a similar decrease for the lensing probability of 
the Ly a spheres. Despite of the baryon and other con- 
tributions we discussed above, we estimate the optical 
depth as 10~ 5 /6 w 2 x 10~ 6 . 

The total number of Ly a spheres across the sky is 
~ 10 11 — 10 13 , depending on the cosmological parameters 
adopted and the redshift range (|Chen fc M iralda- Escudel 
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1200(1 . Since these Ly a spheres have typical diameter 
10 , some of them apparently overlap along the line of 
sight. However, when they are visible in 21cm, they are 
also well separated in redshifts. The density of these 
spheres is between 0.1 to 1 comoving Mpc~ 3 at red - 
shifts 30, c.f. Fig. 5 of IChen k Miralda-Escudl (|200ll . 
while their diameter is about 20 kpc physical, or 600 
kpc comoving, so they are well seperated and can be 
distinguished at different frequencies. As will be shown 
in next section, the optimal bandwidth of observation is 
about 20 kHz, which could easily distinguish these differ- 
ent spheres along the line of sight. There is some chance 
that within the Ly a sphere another first star form, and 
produce another Ly a sphere which overlaps with the 
first, the two spheres would be physically merged. This 
is however not common. The probablity of finding an- 
other halo within a Ly a sphere is about 0.3 (Fig. 8, 
op.cit.), but given the short life time of the first stars, 
the chance of the two stars shine at the same time is 
even smaller. Of course, eventually as more and more of 
these spheres appear, a homogeneous Ly a background 
is built up. At this point, the Ly a flux induced spin 
temperature difference diminishes, and the Ly a sphere 
become invisible. This happens before the physical over- 
lap of the spheres, at a filling factor of about 0.1 (see the 
discussion on the build up of L y a background in §4.1 of 
IChen k Miralda-Escudel (l2006h ) 

Given the probability 10~ 5 to be amplified by a fac- 
tor of 10 or more, there will be at least 10 6 — 10 s of 
Ly a spheres which are strongly lensed in the concor- 
dance model. A radio array with total collecting area 
Atot < 100 km 2 while keeping cost comparable to SKA 
is promising to detect many of them in years survey as 
the discussion bellow. 

4. OBSERVATIONAL FEASIBILITY 

A large [i will significantly reduce the requirement on 
the resolution and r.m.s. noise of a radio survey. The ex- 
act number depends on the detailed array configuration. 
Without loss of generality, we focus on a configuration of 
N x N arrays homogeneously distributed over a square 
with side length L. For an optimal detection, the reso- 
lution 6 p ~ X/L should be similar to the size of the Ly 
a sphere, Here, A = 21(1 + z)cm is the redshifted wave- 
length of the 21cm line. The residual noise per resolution 
pixel (with area 9 2 ) is 
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/cover = Atot/L 2 and Atot is the total collecting 
The bandwidth Av is related to the comoving 
separation r by Av ~ 42kHz[31/(l + z)] 1/2 [d//i" 1 Mpc]. 
For the typical signal ~ — lOOmk and typical comov- 
ing Ly a sphere diameter d ~ 0.6 Mpc (physical di- 
ameter ~ 20 kpc), assuming the integration t — lyear, 
Atot ~ 1 000 km 2 is required to detec t unlensed Ly 
a spheres ()Chen k Miralda-Escudll2006l ) . 

However, a factor of /x amplification in the area relaxes 
the requirement on resolution and thus the requirement 
on L by a factor of /i 1 / 2 . Since the differential bright- 
ness temperature does not change, to obtain the same 



S/N, A to t can be reduced by a factor of (i, if <jt and 
/cover are kept fixed. Thus, for those highly amplified Ly 
a spheres with fj, > 10, a radio array with A to t ^ 100 
km 2 is sufficient. 

A radio array with Atot ~ 100 km 2 can be built 
of fix ed parabolic cylindrical reflectors (Pe terson et alJ 
2006) with reasonable cost. Such design can reduce the 
cost to 5 ~ 10$/m 2 . Given this number, the total cost 
can be controlled to be around one billion dollar, com- 
parable to the estimated total cost of SKA, which has 
unit cost of ~ 1000$/m 2 . This design has large field 
of view (several thousand square degree at z = 30) and 
can scan at least half the sky through the Earth rota- 
tion. Since the field view is large, the telescope does not 
need to point to specific regions of known galaxy clusters, 
although these regions are of primary interest in the pro- 
cess of data analysis. With these advantages, this design 
is appropriate for the detection of Ly a spheres in blank 
sky. 

5. DISCUSSIONS AND SUMMARY 

The cluster which serves as gravitational lens may have 
radio emission, this could complicate the situation. How- 
ever, we expect that the main mechanisms of radiation in 
the clusters-synchrotron and free-free- produce smooth 
frequency spectrum, so we should be able to disentan- 
gle these from the lensed 21cm signal. The cluster radio 
emission does increase the sky background noise, hence 
the boost in the signal to noise ratio is less than ex- 
pected from the lensing magnification. However, radio 
observation of nearby clusters show that not all clusters 
have radio strong emission, if we avoid those with FR- 
II radio sources, then many cluste rs have radio em ission 
power less than 10 25 W/s at 20 cm (|De Youndl2002D . For 
a cluster at cosmological distance (10 3 Mpc), the cor- 
responding brightness temperature of such a cluster is 
much less than the galactic synchrotron foreground, so 
the increase of noise would not be significant. 

Ideally, one wants to de-lens individual Ly a spheres 
to directly measure their intrinsic properties. However, 
clusters with M vir ~ lQ u h~ x M Q at z ~ 1 (Fig. 3) con- 
tribute a non-negligible fraction of the important lenses. 
These cluster are generally sub-critical for sources at 
z r*> 1 — 3, faint in X-rays and SZ, and their mass re- 
construction would therefore be difficult. Furthermore, 
the angular resolution we consider here is comparable to 
the size of the (lensed) spheres. This makes delensing 
very difficult, even if detailed lens information is given. 

In the worst case that delensing of individual Ly 
a spheres is not feasible, their direct detection is still 
valuable. The direct observable, the number of detected 
(strongly lensed) Ly a spheres as a function of redshift, 
is already valuable for understanding the history of first 
star formation and reionization mechanism. By the time 
of successful Ly a sphere detection, galaxy and cluster 
surveys will be significantly advanced. Cluster abun- 
dance and lensing efficiency are promised to be well mea- 
sured through observations of strong and weak lensing of 
galaxies, as well as the X-ray and SZ effect of clusters 
Given these lens information, the prediction of the de- 
tected number of spheres then relies only on the source 
property. Comparing the prediction with observations 

5 Private communication with Ue-Li Pen. 
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will then reveal the nature of the reionization process. 

As a summary, galaxy cluster is a powerfully telescope 
to magnify first star Ly a spheres. Since the size of 
Ly a spheres is comparable to the caustic size of clus- 
ters, the lensed images often merger together. This dis- 
tinctive lensing feature significantly improves the ob- 
servational feasibility of these objects. Clusters with 
M vlr ~ lO 14 - 5 /^ 1 !/© in the redshift range [0.4,1] are 
responsible for most of the strong lensing events (fig. [3]) . 
Regions where these clusters cover should be the primary 
locations looking for Ly a spheres in the reconstructed 
21cm brightness temperature maps. Cylindrical reflector 
array with A to t ~ 100 km 2 and reasonable cost compa- 
rable to that of SKA should be able to detect at least 



10 5 -10 8 strongly lensed Ly a spheres and open a window 
into the very first objects in the Universe. 
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